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ABSTRACT 



An arbitrarily orientable Yagi-Uda array antenna was 
modeled, and a computer simulation run to obtain the input 
impedance, gain pattern and front- to-back ratio of various 
arrays. The model made provisions for the antenna to be 
operated over either a lossy ground plane or aboard a ship 
in seas of specified state. Quick solution turn-around, 
with CRT display, enabled relatively rapid optimization of 
numerous arrays. 

Theory, resultant optimal designs and performances, 
photographs, and program listing are included. 
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TABLE OF SYMBOLS 



Computer 


Description 


ALPH 


Array axis (boresight) elevation 
angle; Input; Variable during 
ship dynamics: see equation (39). 


ALTEM 


A constant; Input. 


K 


Wave number. 


CEE 


Reflection factor; see equation (7) . 


D 


Adjacent element separation; Input. 


DELTA 


Observation elevation angle 
(90-0); Input. 


DLPRI 


P 

Element elevation angle (90-0'); 
Input . 


EPHI 


Electric field; see equation (35) . 


ETHET 


Electric field; see equation (37) . 


EPSLN 


Dielectric constant of earth; 
Input . 


F 


Frequency MHz; Input. 


G 


Gain; see equation (44) . 
G^ is peak. 


VAR 


Sinusoidal angle, 0* -360* . See 
equation (40c) . 


H 


Height of antenna above earth 
plane; Input; a variable during 
ship dynamics; see equation (41). 


HTEMP 


A constant value of height; Input. 


HDBL 


Separation between actual and 
image elements; see equation (33) . 


CUR 


Actual element current. See 
equations (1) , (2) . 


LH 


Element half-length. 



L 


LH 


Element length. 


l p 


LP 


Primary element length. See 
equations (12) . 


L 

s 


LS 


Secondary element length. See 
equations (12) . 


N 


NE 


Number of array elements. 


X 


LMDA 


Wavelength . 


n 


K0S1 


Angle between array boresight 
and observation line; See 
equation (28) . 


* 


PHI, M 


Observation angle; See figure 3; 
Input; M varies from 0 to 90 
degrees. a occurs at peak gain. 

fr 


4>' 


PHIPR 


Element angle; See figure 5; 
Input . 


$ 


PSI 


Angle between element and 
observation line; See equation 
(18) . 


R h 


RH 


Horizontal reflection coefficient; 
See equation (8a) . 


R v 


RV 


Vertical reflection coefficient; 
See equation (8b) . 


V 


RHPRI 


Horizontal reflection coefficient 
used with impedance; See equation 
(8a) , taking 0 = Ot 


V 


RVPRI 


Vertical reflection coefficient 
used with impedance; See equation 
(8b) , taking 0 = 0*’. 


0 


SIGMA 


Conductivity of earth; Input. 


V 


SZ 


See equation (11a) and figure 2. 


V 


SY 


See equation (lib) and figure 2. 


0 


. THET, KAY 


Observation angle; See figure 3; 
Input; KAY varies from 0 to 360 
degrees. 6p occurs at peak gain. 


9 t 


THTEM 


A constant value of 0; Input. 
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e ' 


THEPR 


Element angle; See figures 2,5; 
Input. 




WYE 


Separation between antenna 
reflector and the i^h element; 
See figure 1. 


Y o 


YO 


Y coordinate of secondary 
antenna origin; See equation (10b) 


Z o 


ZO 


Z coordinate of secondary 
antenna origin; See equation (10a) 


z 


Z 


Impedance value. 


[Z] 


ZZPAK 


Impedance matrix; see equations 
(3) and (6) . 



9 



I. INTRODUCTION 



A. BACKGROUND OF THE STUDY 

The most concise equation that describes the radiated 
electric field of an antenna assumes that the antenna 
operates either in free space or over a perfectly conducting 
ground plane. Equations in this category are simple ones 
that have been typically used in text books, and voluminous 
experimental data have been obtained through their use . 
Emphasis has been upon simplification of equations so that 
problems could be solved manually with a minimum of rigor. 
Equations such as these are admittedly inaccurate because 
they omit the component of radiation produced by ground 
reflection or if not omitted the ground surface is considered 
to be a perfect ref lector . ^ 

In reality the simplest antenna operates over a lossy 
ground and established at least a direct and a ground 
reflected wave. When simple free space equations are modi- 
fied for two path propagation over a lossy ground plane and 



l‘The ground wave is composed of a space component and 
a surface component. A particularly important point may be 
made that for the horizontally polarized wave in the VHF/UHF 
range the strength of the space wave greatly exceeds that of 
the surface wave so as to render the surface wave negligible. 
This then allows an accurate field strength to be calculated 
using equations which assume only a two-path propagation 
model--the direct wave and the lossy ground-reflected wave. 
Henry R. Reed and Carl M. Russell, UHF Propagation , (New 
York: John Wiley & Sons, Inc., 1953) , p. 174 . 
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solutions are obtained by computer with a Cathode Ray Tube 
(CRT) display of the gain patterns, interesting results 
accrue. 

B. STATEMENT OF THE PROBLEM 

The initial problem is to obtain and verify the equations 
required to describe the electric field of a multi-element 
linear array, particularly of the Yagi-Uda type. A recent 
work which provided graphic computer solutions to single 
element antennas is the basis of this study, and as such 
this study is a follow-on. ^ Therefore effort centers around 
adapting the existing program to meet the requirements of 
both antenna types. 

C. OBJECTIVES 

The objective of this study is to provide a near real- 
time computer graphic solution of the gain pattern of an 
arbitrarily oriented Yagi-Uda array which is centered above 
two types of planes (with specified t & a) : 

1. the lossy ground plane, and 

2. the lossy ocean plane which rolls and pitches the 
antenna as specified by the sea state. 

Examples of gain patterns and other output parameters are 
illustrated. 



^R. W. Adler and C. B. Robbins, "The Solution and Graphic 
Display of Gain and Patterns for Wire and Linear Antennas in 
the Presence of Lossy Ground", Electrical Engineering Depart- 
ment, Naval Postgraduate School, to be published . 
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D. SCOPE AND LIMITATIONS 

The following assumptions have been made for the study 
of the Yagi-Uda array: 

1. Propagation is confined to two paths. 

2. The ship-ocean model does not make any provision 
to augment the number of wave paths that turbulent 
seas might produce. 

3. The current on the elements is distributed sinusoi- 
dally, e.g. the elements are thin (d <_ X/100) . 

4. The elements may be spaced arbitrarily with arbitrary 
lengths. Thickness can be changed. Assignment of 
the element lengths is constrained by the fact that 
if elements are one wavelength then the solutions 
become indeterminant. 

For purposes of testing the resulting equations, various 
array designs were attempted. The dimensions of the array 
were varied using the method of iterative search to uncover 
optimal horizontal and vertical designs of arrays with and 
without a reflector. 

Optimality is determined by three criteria: 

1. Input Impedance (Z^ n ) 

Where it was possible to do so, the array was 
designed to have a reasonably high resistance, e.g. 

20^ and a reasonably low reactance, e.g. 10ft. 

2. Front to Back Ratio (FBR) : 

Within a satisfactory impedance range, the FBR is 
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3 

maximized. This expression differs in some cases 
from the expression for FBR that is typically used. 

3. Power Gain (G) : 

Finally, within the maximum FBR the gain is 
maximized . ^ 

Photographs are shown of the linear and logarithmic 
results that are obtained from a variety of Yagi-Uda 
designs placed at different heights above ground and ocean 
environments . 

Preliminary tests were made which showed that gain 
varied with change in element thickness, however for the 
results obtained throughout this study the element thick- 
ness to element length ratio remained fixed at 1/200. 

Detailed operating procedures are found in Appendix A. 



^For purposes of the study the following expression 
for FBR isused 



FBR = G(9,4>)max ® 1. 4> ^_180 

G( y ,<t)max , 1801 4 1360' * 



This differs completely in some cases from the expression 
that is typically used: 



FBR 



G(6 ,90) 

gt rrrm ' 



^The power gain used is measured with respect to an 
isotropic source. 
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II. THEORY OF THE MODEL 



A. ARRAY IMPEDANCE AND ELEMENT CURRENT 

The radiation problem requires the solution of the 
individual element currents, which will become the basis 
for calculating the gain of the array. The assumptions 
are that the driven element is excited by one volt, and 
that the antenna is thin and therefore the element currents 
are sinusoidal. The equations for the current column 
matrix and for I ma x i are 

[i] = [z]" 1 [V] (1) 



L max x 



sin g l i 



( 2 ) 



[z] is the combination of the free space impedance and 
image impedance which is found by 



[Z] = [Z] + C [Z * ] 



( 3 ) 



[z] and [Z 1 ] are NxN dimensional — [Z'] being the mutual 
impedance matrix relating the actual and image elements as 
indicated in figure 1. It should be noted that the matrices 
are complex. Complex inversion presents no particular 
problem except that the [z] matrix must be arrayed according 
to the basic arithmetic operations that follow: 
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Figure I. Side view of a 3 element 
Yagi-Uda array tilted 
degrees abovethe ground. 
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[V] 



(4) 



([R] + j [X] ) • ( [ I R ] + j C I x ] ) = 



which can be written as two equations 



[R][Ir] 


- [X][I X ] 


■ Mrs 


[X][I R ] 


+ 

1 — 1 

& 

1 1 

1 — 1 
H 
H 

i I 


■ to*. 



(5) 



which can be solved according to equation (1) , where 
equation (3) now appears as: 



[Z] 



[r] 


[-X] 


[X] 


[R] 



(6) 



[z] is a 2 Nx 2N matrix. An array of maximum size (5 

elements) requires inversion time on the SDS-9300 computer 

of 5 seconds--roughly one-eighth the total solution time 

needed for the complete aximuth and elevation patterns. 

The reflection factor C in equation (3) is a function 

of the angle which the antenna element makes with the 

5 

ground. C is written as 



C = e 3 A '(R h ' CO s A' + jR v 'sin A'). 



(7) 



The reflection coefficients R^' and R v ' are the values 
obtained when 0=0* is substituted into equation (8) . Theta 



5 

M. T. Ma and L. C. Walters, Power Gains for Antennas 
over Lossy Plane Ground , ESSA Technical Report. ERL 104-ITS 
74 (U. S. Government Printing Office, Washington, D.C., 

1969) . 
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is equal to zero because the coupling which takes place 
between the free-space and the image element (thereby 
producing a mutual impedance Z 1 ) occurs with the image 
directly beneath the actual element. Equations (3) and (7) 
show that the horizontally polarized field (A'=0°) and the 
vertically polarized field (A'=90°) give weight to the 
value of [ Z ' ] by the values C=R^ ' and C=R v ' respectively. 
Since the solution to equation (3) is actually independent 
of 0 and <f> , one solution for this equation satisfies the 
gain expression at any position of observation. 

The values of R^ and R v , which are also used in the 
gain equation, are a function of the observation angle 6 
as follows:*’ 



cos 0 



(8a) 



*v 



COS 0 



COS 0 



COS 0 



+ 



+ 




B 





(8b) 



where 



A 



22 

or 

‘ 22 



[1 - ( 



. . 2 n 2 
sin e) ] 

b 2 



2 [ e r 

B[e r 



j 

j 



U)G 



s-i 



o 



1.8 a 10 



'MH 



% 

] 



(9a) 

(9b) 

(9c) 



6 Ibid . 
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The matrix [z ' ] is identical to equation (6) except 
that it is written in terms of [R'] and [X']. The geometric 
orientation of the antenna elements are central to the 
solution for the mutual impedance matrix [ Z ' ] . Figure 2 
shows an arbitrarily oriented single element and its 
corresponding image in the ground plane. Although the 
true coordinates of the antenna are in the xyz coordinate 
system, the equation for solving for the impedance fixes 
the primary element along a vertical axis — here it is shown 
to be z''. The secondary element becomes the image element 
and is separated from the primary element by the distances 
Z Q and Y . The linkage between the two coordinate systems 
xyz and z'y'z' is 9'. The equations for Z Q and Y Q are 
observed to be 



II 

0 

[SJ 


-2h cos 9 ' 


= -2h sin A' 


(10a) 


II 

o 

>H 


2h sin 9 ' 


= 2h cos A' 


(10b) 


S 2’ - 


S cos 29 ' 


= S cos 2 a' 


(11a) 


V = 


-S sin 29 ' 


= -S sin 2 a' 


(lib) 



The solution of self impedance is performed by assuming 
that Z q =0 and Y 0 = 2 r, where the radius r is assumed equal 
to i^/200. (This makes £^/d = 100, which for SL^ = A/4 
satisfies the thin antenna specifications of d < A/100) . 

The mutual impedances of the free-space array, where the 
elements are all parallel, are solved using Z Q =0 and Y Q =y^, 
the separation between the two elements in question. Under 
dynamic ship motion, which affects 0' and h, the antenna 
height above the earth plane, the mutual impedance 
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Figure 2. The geometry of an arbitrarily oriented element 
in free space and its image (used for solution 
of mutual impedancesZ') 
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matrix [z'] will be affected according to equations (10) and 
(11) , but the self and mutual impedances of the [z] matrix 
will not be affected since it represents the free-space 
array impedance . 

Impedances are solved according to the following:' 



l s/2 

R = : 3 o / jf 

-V 2 



s “ (2 “ rl) - (Z 0+ S 2 . + %, + (Z^-^) 



r l 



2T 



2X 



cos (£g) (Z 0 +S z ')1 
2X J 



c • 



Y+Sy' 



2 sin(2TTr) cosf^p) - sin(2^r]_) _ sin(2Trr 2 ) 
r X r i r 2 



• S„ 



sin (2 t£e - | s | 
2X 



ds 



X = -30 



L s /2 ( 

JJ 

cos ( 2 tr 2 ) 



cos( 2 ' ,r i ) <z„+s 2 '+^£) 

r l 



(Z 0 +S 2 *-^) - 



o‘ r °2 ^Jx' 

2 cos(2iTr) 



r 2 

(Z 0 +S z ' ) 



ttL_ 
cos ( E) 



2X 



y + 
' 1 Z 



(Y 0 +Sy') 



-cos(2ur^) cos(27rr2) 

r l ^2 



2 cos(2TTr) .ttL,-,. 

COS ( P) 



L s 



^ s in ( 2 tt <2 ^ I s | ) 



(12a) 



ds 

(12b) 



^H. C. Baker, and A. H. Lagrone ; "Digital Computation 
of the Mutual Impedance Between Thin Dipoles;" IRE 
Transactions on Antennas and Propagation; March, 1962; 
AP-lO , No. 2; pps 172-178. 
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V 







[p 2 + (Z 0 + S z '-iE) 2 ] 



(13a) 



(13b) 



(13c) 



r 



2 



(13d) 



Where the distances r, r^ , and r^ , are the respective 



distances in figure 2 from the bottom, center and top of 
the primary antenna to the differential, ds, on the 
secondary antenna. 

B. THE FAR FIELD 

To solve for the electric field, a point of observation 
is chosen, and each element in the array is viewed. Each 
element contributes to the total field according to its 
vector phase and amplitude, its separation from some 
reference element, the tilt angle of the array, and the 
bearing to the point of observation. Figure 3 shows the 
geometry of a 2 element horizontal Yagi-Uda tilted a degrees. 

The equation of the field produced by one element located 
at the origin is first developed. Then, the presence of 
other elements in the array is taken into account. 

The equation at the point P for an element at the origin 



Current flowing in the incremental length dx produces a 
field, and there will be a phase associated with it 



is 





sin[6(£ i ±x)] e 



(14) 
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z 




Figure 3. Yagi-Uda array with horizontal elements and 
a tilted axis, as viewed from a distant point P. 

(<£'=0,0=90) 



with respect to point P as the differential length moves 

from to + . The delay is x cos \p , which is found 

2 ^ 

by the dot product of OX and OR, as follows: 

OA = OX • OR 



where. 



II 

<o 


x [i sin 6 ' 


A A 

cos 4> ' + j sin 6 ' sin <p ' + k cos 0 '] 


( 15a) 


* 

OR = 


A 

[i sin 6 


A. 

cos + j sin 6 sin $ + k cos 0 ] 


( 15b) 


so that 








OA = 


X [sin 6 


sin 0' cos + cos 0 cos 0'] 


(16) 


or 








OA = 


X cos l p 




(17) 


where 








cos iJj 


= sin 


0 sin 0' cos ( 4> — <f> 1 ) + cos 0 cos 0' • 


(18) 


It is 


now clear that s = r - x cos and 




I . = 

1 


I sin 

o . 

l 


r n f r- X COS ♦ )1 

c B ( £ ± ±X) ] e J Lt ( C )J 


(19) 


where 

£±X = 


M + x 
\ a - x 


X X 

V A 

O O 

% 


(20) 



Although there is a phase difference between r and s, at 
points taken far from the origin their magnitudes are 
approximately equal. 

Observe that by translating the xyz origin as well as 
the off-axis element to the point B located along the 
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observation vector r, the magnitude of the new observation 
vector relative to the new origin is BR, which differs in 
length from the old observation vector OR. The amount of 
change equals that which affects the observation vector s. 
Furthermore the angle by definition remains unchanged. 
Therefore, the expression r - s = x cos 4> holds for every 
radiating element in the array, not just for the element at 
the origin. From this it can be concluded that the integral 
equation for the field is identical for each element in the 
array. 

O 

Following the method used by Kraus: 



£ . 



j 60 tt sin <t> 




I . dx 



1 



( 21 ) 



. . 2 
Substituting for I and letting B =. — 



c 



we see that: 



E<{,. = j— 3P — ^ sin <p I 

1 T 



r 



o. 

i 




sin(B*U±Bx)e 



jBxcos* ^, 22 , 



When this expression is integrated the result is: 




[cos(6£^cos ^)-cos(B£^) ] . (23) 



• 2 , 

r sin \p 



®J.D. Kraus; Antennas; (New York: McGraw-Hill Co. 1950) 
pp. 135-141. 
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For an array of elements the field is the superposition of 
the fields of each element taken separately. is then 

the sum of E<j>^ given by 



N 

_ e~3 r 60 sin ^i \ [cos (8£-cosiM -cos (0£^) ] . (24) 



r sin 



~ / j “i 

i=l 



It is recognized that the time term, e~^ cot has been dropped 
and that 8 has been substituted for w /c . 

A 

The array factor, or the phase delay OB introduced by 
the separation between elements, is next taken into account. 
From figure 3, using the reflector as reference, the 
separation between elements as viewed from point P is 



OB = OY-OR 


, where 










A 


* 


A 








OY = 


jyicos 


a + ky^sin 


a 




(25) 


A 


A 


A 


A 






and OR = 


i sin 9 


cos<j) + j sin 


6sin<j> + k 


cos 6 


(26) 


gives 


OB 


= y^(cosa 


sin 0 sin 


ip + sin a cos 0 ) 


. (27) 


If cos 0 


= cos 


a sin 6 sin 


. 4> + sin a 


cos 0 


(28) 


then 


OB 


= y^ cos 0 




(29) 


represents 


the phase lead of 


the i^h element. 





Having found the slant separation between elements to 
be expressed by equation (29) , the current for any element 
in the array may be written as 

I = i e j8YiCOs n . (30) 

Oj_ Max i 

The phasor sum of currents in the array becomes 
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(31) 



I 



o 




jBy^cos U 

I • e 
max 1 



Upon substitution of equation (31) into equation (24) the 
result is 




2 



, n -jBt . , 

60 e sin <p 

. 2 

r sin $ 



E 



-r j0y-:cos 

• e 
max l 



[cos(B^ *cos 'M - cos (6®-^) ] (32) 



where I max ^ is found from equation (2) . 

Equation (32) represents the free space radiation of 
the horizontally polarized Yagi-Uda. The equation does not 
consider the ground reflections. Figure 4 shows the direct 
and ground reflected waves, and the geometry shows the dif- 
ference between their path lengths. The element separation, 
written as: 



2h^ = 2 (h + y^ sin a) (33) 

A 

is a function of the path difference vector AB, since 
AB = 2h^cos 0 . 

Thus 

E <J>T = E <f> + R h E <P e ~ j2 * hiCOS 9 . (34) 



And finally, equations (32) and (34) are combined to produce 
the final space wave in the phi direction: 



Ee 



j 60 e sin 

• 2 . 
r sin i|> 



max; 



[1 + R h e - j28h i cos V 

e 88 Yi C ° S r t COS (6£ iCQSi^) - cos(g£^)]. (35) 
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Figure 4. The tilted Yagi-Uda receiving horizontally 

polarized waves via direct and reflected paths. 

(The array is tilted to present less confusion 
in distinguishing angles) 
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Observe that when cos fi = sin 6 sin <j> and h^=h equation 

(35) agrees with Ma and Walters equation for the horizontal 
9 

Yagi-Uda. The equation for the 6 component of the horizon- 
tal Yagi-Uda, E 0/ is not derived with the same rigor since 
it represents the cross polarized radiation, which exists 
only off-axis and is of minor interest. The equation is: 



_ i 60 e J cos 0 cos 6 

= 71 

r sin <|) 



r — j 2 Sh-: cos e 

> I [1 - R e J 1 ] 

/ ^ max- ” 

i=l ± 



v 



e j B Yi cos fi [ cos(6 £ iCOS <p) _ costBAi)] . (36) 

10 

The equation for the vertical dipole is well known. 

If R v is used insteady of R^ then equations (31) and (34) 
may be applied to the standard free space equation to arrive 
at the general solution for the vertical array: 



' 9t 



j 60 e 



-j8r 



N 



r sin 6 



I [l + r e -j2Bh iC os 6-, 
i max: L1 + K v e J 



i= 



DSYiCO ^[ CO s(8A^cos '!’) - cos(8&^)] (37) 



Specifically, cos <P = sin 6 sin <P and cos ft = cos 6 
since a = 6 = o for the vertical Yagi-Uda. 



o 

M.T. Ma and L.C. Walters, op . cit . p. 41. 
l^Kraus, loc.cit. 
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C. THE SHIP-OCEAN MODEL 



Figure 5 shows the effects of the dynamic ship model 
upon the Yagi-Uda. Changes in the sinusoidal waves, 
considering that a roll of eight degrees per sea-state is 
produced and a bow pitch of three degrees per sea-state is 
produced (representative of a light cruiser) , will affect 



the angles e ' and a as follows: 

0' = 6 1 ' ~ (A^cos <j> * + A 2 sin <j> ' ) (38) 

a = a t - A^sin <j> ' + i 2 cos <j> ' (39) 

where 

A^ = (wave) sin (course) (40a) 

A 2 = (wave). 3 cos (course) (40b) 

and where 

wave = 8 (sea state) sin (y) . (40c) 



The effective height of the array will change according 
to 

h = h fc cos A^ cos (41) 

It should be noticed that the parameter 4> ' until now 
has been somewhat extraneous for the solution to antennas 
over land, since variation of the observation angle <j> ' and 
the angle <j> accomplish the same thing. However the angle 
<f> ' now is significant in specifying the antenna orientation 
with respect to the ship heading, which then determines the 
aspect that the wave presents to the antenna. 



29 




Ships 

Heading 



Figure 5. Orientation of the Yagi-Uda 

aboard a ship. Angles a, and a z 
represent roll and bow pitch 
respectively. 
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D. GAIN 



The gain of 


an antenna is equal to 


the ratio of 


the 


power intensity 


to the 


power density. 


and is expres: 


sed as 


follows : 


G 


_ 4$tW r 

w. 

m 




(42) 


where 


w r 


= r 2 ( E x H*) 


watts/solid 


angle 


or 


W r 


2 2 
= r |E| 

120* 




(43a) 


and 


W. 

in 


= I I- I 2 R. 

1 bi in 


watts . 


(43b) 



The final expression for gain appears as 

G M r 2 |Ee T | 2 + 1 E 4 > t | 2 , 

30 I Ibl 2 R in 

The equation for gain becomes independent of r when the 
equations (35) and (36) are substituted into the equation 
above . 
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III. RESULTS OF ARRAY PERFORMANCE 



A. RESULTS OVERVIEW 

Basically, four types of arrays were examined for 
optimization: horizontal and vertical with and without 

reflectors. Experiments were conducted to determine the 
design of arrays that afford optimal performance as defined 
in the Scope and Limitations section of the study. Starting 
with the two element array, the parasitic element function- 
ing both as a director and as a reflector, the element 
spacings and lengths were varied to obtain optimization. 

The arrays were lengthened by adding one element at a time, 
up to a total of five elements, each time solving for an 
optimal design. Numerous iterations were required, and the 
time required to optimize an array varied from ten minutes 
for a two element array to slightly over three hours for a 
five element array. (The vertical array was more difficult 
than the horizontal, and the array without reflector was 
more difficult than that with reflector.) Not all of the 
multi-element arrays were optimized according to the rigid 
criteria chosen for this study. Specifically, the vertical 
designs having four and five elements were arrived at by 
relaxation of the impedance threshold criterion thereby 
producing a suboptimal design. 

B. THE TWO ELEMENT ARRAY 

Referring to Table I for the horizontal two-element 
array with a director, maximum FBR is maintained when, for 
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TABLE I 



OPTIMAL 2 -ELEMENT HORIZONTAL ARRAY 
AT VARIOUS ELEMENT SPACINGS USING 
THE PARASITE AS A REFLECTOR AND AS 
A DIRECTOR 



(Gain and Front-to-Back ratio in db) 
(Spacing and length in X) 




f =146. h=lm. In both cases the length of the driven 
MH 

element is X/2. (L^ and L are director and reflector 

lengths respectively.) (The subscript p denotes peak, 
e.g. the length of the parasitic element was varied until 
the peak FBR was first obtained, then within that value, 
the gain peak was found by making finer adjustments to 
length . ) 
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increasing element spacing, there is a corresponding decrease 
in director length. Apparently, best performance occurs when 
spacing is .01X (,021m) because this is where the maximum 

FBR appears. However, no mention has been made about the 
affect upon impedance of varying the element spacing. 

In general it was found that the resistance and inductive 
reactance decreased together (from 70 to 40 ohms and from 
40 to 10 ohms respectively) as elements were placed closer 
together. When the spacing closed to . 05X, the currents in 
the elements evidently became very large, and the value of 
resistance became less than one ohm (too low for convenient 
matching) . Thus the FBR and G are insufficient criteria of 
optimality since impedance can fall (in some cases may rise) 
to unreasonable values. Optimality in the case of the di- 
rector array in Table I actually occurs at .IX (.205m) 
spacing . 

When the parasitic element is used as a reflector it is 
not necessary to change the reflector element length as the 
spacing between elements is changed in search for the 
maximum FBR. The maximum FBR occurs at . 5X (1.03m) for all 
spacings except at .25X (.515m); that is ,the maximum FBR 

occurred when both element lengths were the same. 

The gain curves shown in figure 6 differ from other 
published experimental and theoretical data (and there are 
fairly wide variations in these results too) , inasmuch as 
the value of the peak gain is somewhat higher than that 
typically recorded. This is because published data typically 
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Power Gain over Isotropic db 




Figure 6. Comparison of gain vs 
element spacing for optimal 2— 
element horizontal arrays. These 
curves represent data in table I , 
where the array is a reflector- 
driven element combination inone 
case and a director-drivenelement 
combination inthe other case. 
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does not reflect a lossy ground plane (two-path propagation 
is ignored). Also, the gain of the reflector array exceeds 
that of the director array, which generally disagrees with 
other findings. Again, the explanation for these differences 
is tracable to the assumption that a lossy ground wave is 
being propagated along with the direct wave.^ 

C. THE MULTI-ELEMENT VERTICAL AND HORIZONTAL ARRAYS 

Tables II and III list the optimal and suboptimal design 

parameters for the f m ^=30 arrays. Notice that the four and 

five element vertical arrays have relatively wide reflector 

spacings — these are suboptimal designs. By comparison with 

the optimal design that was found for the four element 

vertical array (where R^ n >_ 20ft) the FBR is 16.4 db whereas 

the suboptimal design (Rj_n < 20ft) produced a FBR of 22.4 db 

12 

— quite a drop. The two designs differed completely. 
Furthermore they responded to slight parameter variations 
differently . 

The four element vertical array was fine-tuned for 
maximum FBR only to discover that the impedance was 
unacceptable; all subsequent fine adjustments, introduced 
for purposes of raising the resistance to a value above 20ft, 



R. Fry, The Yagi-Uda Aerial — A Short Design Review 
and Bibliography , ( Valcartier , Quebec : Canadian Armament 
Research and Development Establishment, May 1966), p. 23. 

■^The optimal design, as compared with the suboptimal 
design shown in Table II, was Lj_= . 4 9 , . 47 , . 4 5 , . 4 3 ; d^=.22, 
.19,. 15 (L and d in X). 
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THESE DESIGNS ARE ACTUALLY SUB0PTIMAL 



RESULTS OF OPTIMAL 30 MHz HORIZONTAL AND VERTICAL ARRAYS WITH REFLECTOR AND OPTIMAL 
HORIZONTAL ARRAY ORIENTED VERTICALLY 
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caused a drastic reduction in FBR — below 16 db. Thus, this 
particular design was abandoned; another completely different 
design provided an acceptable impedance with a peak fine- 

tuned FBR value of 16.4 db (compare parameters in footnote 

/ 

#12 with the comparable parameters in Table II) . 

Except for the two larger vertical arrays (four and five 
elements) , all changes that were made for the purpose of 
increasing Rj_ n or reducing X^ n inevitably improved the FBR 
figure and in none of these cases did the manipulation of 
parameters for purposes of adhering to the impedance 
criterion result in obtaining a lower FBR or G value. 

Because such manipulation did, however, cause a drop in the 
FBR for the two vertical arrays they are considered 
anomalous . 

13 

The optimal horizontal and vertical arrays have 
tapered elements in the region next to the driven element, 
but the last few end-directors of the four and five element 
horizontal array show a slight inverse taper. Writers, 
however, generally agree that elements and spacings that are 
gradually tapered towards the end of an array will usually 
give best results. 

It is obvious that the optimum parameter measurements 
obtained for the reflector-driver combination (two-element 
array) do not necessarily ensure optimization of successive 
experiments involving the addition of another director. A 



-*-^The optimal arrays include all horizontal plus the two 
and three element vertical arrays. 
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design readjustment of all parameters is required each time 

the array length is changed (e.g. each time an element is 

added or removed) . A corollary to this for arrays with more 

than two elements is that the FBR will decrease with the 

addition of a director whose length and spacing are identical 

to those of the preceding director. This was observed in 

14 

testing the uniform array at f MH =10 where the FBR dropped 
from 14.6 db as a four element array to 10 db as a five 
element array. 

With regard to director spacings, perhaps it is more 
than just coincidental that the designs which produced 
optimization of the longer horizontal arrays resulted in 
equispaced director elements. Nothing could be obtained 
to substantiate whether other analysts would agree. 

From Table III, comparison between the horizontal three 
and four element arrays as well as comparison between the 
four and five element arrays shows that the FBR improvement 
is greater than the G improvement as array lengths increase. 
An added advantage, besides increasing FBR and G, is that 
the greater the number of elements used in the optimal array 
the greater the tendency for automatic reactance cancella- 
tion. While the reactive components of the two and three 
element arrays are small they did not become so as a 
coincidence of optimization; manipulation of parameters was 
necessary to obtain reactance decreases in conjunction with 



14 Design: L i = . 533 , . 500 , . 434 , . 434 , . 434 and d ± = .233 , .100, 

.100, @ h=3 (L, d and h are in A) . 
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the need to increase resistance (and, as inferred earlier, 
these manipulations in the case of the two and three 
element arrays produced slight decreases in FBR and G^) . 

Comparison between the optimal horizontal array and its 
vertical counterpart indicates that the horizontal array 
consistently out-performs the vertical array in FBR, but 
matches it in G. Furthermore, when the optimally designed 
horizontal array is operated in the vertical position and 
vice versa a degradation in performance results. An improve- 
ment of almost 16 db in FBR and 1.5 db in G can be obtained 
when the array is optimized in the position for which it is 
intended to be operated. 

Table IV gives the results of operating the optimal 
array without its reflector; it is interesting to observe 
the importance that the reflector has in determining the 
overall array performance. Practically no directivity was 
obtained without a reflector even with a four element 
(three director) array. The unusually large G that results 
for the three and four element vertical arrays are anomalies 



l^Fishenden reiterates the popular belief that a 
convenient method of altering the input impedance without 
affecting FBR or G is by varying the spacing between the 
reflector and the driven element from A/8 to X/4. Such 
attempts in this study did however noticably affect the 
FBR and G (particularly in the case of the vertical four 
and five element arrays, which v/ere particularly trouble- 
some) . R. M. Fishenden and E. R. Wiblin, "Design of Yagi 
Aerials", Proceedings IEE, 56, Pt. 3, (Jan 1949), p. 6. 
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TABLE IV 



THIRTY ARRAYS DESIGNED AS SHOWN IN TABLE II 
BUT WITH REFLECTORS REMOVED 



1 

Horiz @ Horiz i 

Orientation 


i Vert @ Vert 

J Orientation 


N 1 


2 


3 


4 


1 


2 


3 


4 


1 

1 

G | 
t 


7.8 


11.5 


5.0 


7.9 


7.8 


4.8 


15.5 


15.6 


1 

FBR | 

1 


0 


0 


1.7 


1.8 


0 


0 


2.8 


2.9 


R in' 

1 


60.6 


29.4 


4.1 


33.0 


60.6 


34.2 


18.1 


16.1 


x in' 

1 


-12.3 


-0.1 


-15.3 


-9.0 


-12.3 


-9.7 


-6.4 


-4.2 
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and cannot be explained (undoubtedly it is connected with 
the fact that these two arrays are actually subo’ptimally 
designed to begin with.) 

Such poor performance is not without remedy, however, 
as the results show in Tables V and VI. Arrays without 
reflectors may also be optimized with considerable improve- 
ment in performance. Unlike the arrays with reflectors, 
the optimal design of arrays without reflectors gives 
identical optimal results regardless of whether the array 
is operated in the horizontal or vertical position. Inas- 
much as the reactance is much greater than 10ft, these 
arrays do not represent optimal designs according to the 
purposes of this study, yet none better could be found. 

D. SELECTED INTERESTING PATTERNS 

Numerous photographs were taken of the following antenna 
patterns (in a few cases comparisons between the linear and 
the log gain plots are shown) : 

1. a three element array operating over dry land, fresh 
water and sea water, 

2. a three element array at different heights over land, 

3. several three element arrays operating on a ship 
underway in heavy seas, 

4. a three element array with its axis tilted at various 
angles from the horizon and 

5. some resulting patterns when an array is operated 
off-frequency . 
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TABLE V 



PARAMETERS OF OPTIMAL 30 MHz ARRAYS WITHOUT REFLECTOR 

(Horizontal & Vertical Optimal Arrays Identical) 

(L and d in A) 




TABLE VI 

RESULTS OF OPTIMAL ARRAYS WITHOUT REFLECTOR 

(Specifications and parameters same as in Table V) 
(G and FBR in db, R and X in ohms) 





N 2 


3 


4 


G 


19.7 


13.6 


10.1 


FBR | 


0 


14.9 


26.2 


R 1 


4.8 


76.6 


73.9 


x ! 


-6.2 


57.4 


74.5 
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Figure 7 Optimal 30 MHz 

3-element horizontal array: 
linear plot (a) and log 
plot (b) , 16 

G p = 12.6 db 
FBR p = 22.0 db 
Rin = 31.6 Ji X in = .1JL 




Figure 8 Optimal 30 MHz 

4-element horizontal array: 
linear plot (a) and log 
plot (b) . 16 

G p = 13.3 db 
FBRp = 23.9 db 
R in = 36.1 a X in = -2.3-fl- 




Figure 9 Optimal 30 MHz 

5-element horizontal array: 
linear plots la) & (b) , and 
log plot (c) . 6 

Gp = 13.5 db 
FBRp = 27.3 db 
R in = 28.3^. X in = -1.0-n. 



16 

Parameters that are not shown in Table II but which . 
are common to all the 30 MHz arrays are: ht=3, e=5. a=10 - 

<{> =90. 6 =85, however figure 9 (c) differs in observation, 

with <(> = 20? and 6=60. 
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(a) 



Dry 

f!r p 



Land 

= 10.7 db 
= 9 . 6 db 




(b) 

Fresh Water 
G p = 9.9 db 
FBRp = 9.5 db 




(c) 

Sea Water 
G = 12.7 db 
FBRp = 9.5 db 



Figure 10. Ten MH three element vertical 
array radiating over three types of terrain. ' 




Linear 




Figure 11. 150 MH Optimal 

three element horizontal 
array — linear plot (a) , and 
log plot (b) . 



pippi 

mmm 



fIr. 

Rin 

^in 



* 12.9 db 
= 22.0 db 
= 28.7 Jl 
=-13.5 jt 









(b) 

Log 



a=10 

e=80 



17 Li =.506, 

^ r <f ) p = 90. 

. , a = l0”2 f 



18 Li=.510, 
4>p = 9 0 . ; 0 p =89. 



500 


,.4 


50; d i= 


.233,. 133; h t =2.; 


f MH 


=10. ;e= 


10. , 


0 p 


= 84 


. (land) 


; £—80 . , o 2 . , i) |> 


=90. 


, 0p=7 6 


. (H 2 0) 


*p~ 


90. 


, 0p=75 


. (sea) . (L, d anc 


1 h t 


are in 


A) 


490 


/ • 4 


60; di= 


.250,. 08; h t =15; 


e=5 . 


; a=l 0 - 


3. 

/ 


(L, 


d 


and h t 


are in A) 
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ic; 

h t = 1.0 

G p = 12.9 @ p =76 . 0 

FBRp =27.0 



(d) 

h t =2.0 

G = 13.3 @ =83. 

FBRp =27.3 v 



Figure 12. Ten MHz 3-element vertical arrays (a) & (b) 
and horizontal arrays (c) & (d) at different height 
over land. 19 (g and FBR in db) 



l^see figure 16 for parameters. (h t is in A.) 
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(a) 

a) = 0 . 

stable 




(b) 

a) = 30 . 

pt.roll = 15.4° 
up pitch = 5.5® 

0' = 71.5 
a = 64.0 
h = 1.82 




(c) 

u = 90 . 

st. roll = 31.0* 
dn. pitch = 11.0° 
6 * = 120.0 
a = -11.0 
h = 1.68 



Figure 13. Ten MH three element horizontal array 
aboard ship in rough seas. Sea state = 7 and relative 
^Vopi-ion of waves = 40. Gain varies. 20 




u = 0 . 
stable 




vdj 

id = 30 . 

pt.roll = 21.0° 
up pitch = 3.6° 

8 ’ = 69.0 
a = 3.6 

h = 2.80 




(c) 

w = 90 . 

pt.roll = 41.5° 
up. pitch = 7.2* 

0' = 48.5 
a = 7.2 

h = 2.20 



Figure 14. Thirty MH three element horizontal array 
aboard ship in rough seas. Sea state = 6 and relative 
direction of waves = 60. Gain varies. 2 ^ 



20 Li=.506, .500, .450; di= . 233 , . 133 ; e=80.; o= 5.; $=90.; 
0 =80.; h t =2. (L, d and ht are in X) 

21l • = . 51 , . 50 , . 46 ; di=.25,.08; e = 80.; o= 5 . ; <j>=90.; 
0=85.; h t =3. (L, d and ht are in X) 
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= 0 . 
stable 




15 . 4* p. roll 
5 . 5°u. pitch 
0' = 58.5 
a = 5.5 

h =14.5 



31 . 'p. roll 
11 u .pitch 
9' = 59.0 
a 11.0 
h =12.7 




= 240. 

27 .* s.roll 
9 . 6°d. pitch 
0 ' = 120 . 
a = -9.6 

h = 13.2 



Figure 15. 150 MHz 3-element horizontal array aboard 

ship in rough seas . Sea state = 6 , and relative direction 
of waves = 40. Gain varies. 22 



^Other parameters include L^=.51, .50, .46; d^=.25,.08; 
"£=80., a=5., <f>=90 . , 6=89., h=15 . 
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o o 





a 



P 



90.0 

3 . 506=76 . 0 
7.7§e =1.0 




(e) 

a =90.0 
G =13.3 
FBR =0.0 
♦ = 180.0 

0 = 76.0 




(f) 

a =90.0 
G =13.3 
FBR =0.0 
<P =0.0 
e = 50 r 



Figure 16. Ten MHz 3-element horizontal array at 
various tilt angles of array axis (bore sight). 23 



d 

h 




^Various antenna parameters are: L^=.506 

233, .133; e=10., a=10“3, = 9 0 . , 0=76. for 



, .500, .450; 
(a) - (d) , 



50 









Figure 17. Ten MHz 3-element 
horizontal array operating over 
land at f M H=30. (For antenna 
specifications see the ship-ocean 
model figure 14.) 




Figure 18. Thirty MHz 3-element 
horizontal array operating over land, 
above and below frequency for which 
array was designed. (For specifications, 
see the ship-ocean model figure 14.) 



51 




Figure 19. Ten MHz 5-element 
horizontal array over land. 
Comparison of linear (inter) 
plot with log plot. 

G p =14 . 6 db, FBR p =6 . 5 db. 24 




Figure 20. Ten MHz 3-element 
horizontal array over land. 
Comparison of linear (inter) 
plot with log plot. __ 

G n =12 . 2 db, FBR =17.7 db. 

p p 



24various parameters include: L-=.534. .500, .466, .466, 

.466? di=.25, .167, .167, .167? o =10. , o =10 -2 ? <f> =90. ? 6 =77. ? 
h=l . 

^Parameters for this array are: L^=. 550 , . 500 , . 433 ? 

d^=. 233 , . 133 ? all other parameters are the same as in 
footnote 24. 
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E . SUMMARY 

There is a somewhat higher gain and an increasingly 
higher f ront-to-back ratio as the element spacing widens 
in the case of the optimal two element horizontal array 
possessing a reflecting parasite as compared to the optimal 
array with a director parasite. 

The vertical multi-element array will produce a high 
f ront-to-back ratio — with values in the neighborhood of 
those of the horizontal array; however, if resistance is 
to be maintained above 20ft it appears that the resulting 
optimal vertical array will not provide as high a front-to- 
back ratio as the optimal horizontal array. The optimal 
design of a multi-element horizontal array provides both 
a maximum f ront-to-back ratio and a near resonant input 
impedance, whereas the best design of a multi-element 
vertical array represents a compromise between these factors 
e.g. a trade-off sacrifice exists. 

The optimal design arrived at for a horizontal array 
becomes a suboptimal design when operated vertically, and 
vice versa. 

Upon adding or eliminating directors from an array which 
has been designed for optimal performance — with the stipula- 
tion that parameters of all directors be kept uniform, the 
performance of the array will no longer be optimal. 

Operation of the horizontal array ©ver sea water environ 
ment provides a higher gain with respect to operation over 
land, but over fresh water provides a lower gain than over 
land . 
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Gain is considerably higher in the case of the vertical 
array and slightly higher with the horizontal array when 
they are radiating from a height of 2A as compared with IX. 

Gain varies considerably when the platform upon which 
an antenna operates is in heavy seas. 

The gain of a horizontal array aimed vertically into 
the sky is lower than when it operates with no tilt angle 
whatsoever . 
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IV. RECOMMENDATIONS 

At the beginning of this study the existing program 
which had been written for eight different single element 
antennas permitted the compilation of, at maximum, only 
three different antennas — the computer memory was exceeded 
when more options were added. Now the revised program 
allows all nine antennas to be compiled, with the eight 
unused antenna subroutines residing externally on drum. 

This permits consideration of further expanding the main 
program capabilities inasmuch as the effective computer 
memory has more capacity. 

In regard to expanding the program's capability for 
solving array problems, it is recommended that the graphic 
input portion of the program be extended to allow for more 
element length and space parameters. This would permit 
study of long arrays (arrays having a cumulative spacing 
that total more than one wavelength) . 

It is recommended that an optimization routine be added 
to the main program which would encompass one or more of 
the following criteria: gain, front-to-back ratio, 

resonance of the driven element (Z 22 ) • thickness of the 
element, and impedance (or standing wave ratio) . For this 
purpose, a more powerful computer — which gives faster 
solution time — might be preferable. 
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APPENDIX A 



DESCRIPTION OF PROGRAM OPERATING PROCEDURES 

The program for the system of antennas as originally 
written has been changed to accommodate the Yagi-Uda. 
Accordingly, the program was rewritten in part to obtain 
more effective use of limited information storage in the 
computer. Redundancies between antenna programs have been 
reduced, the shipboard simulator and the gain pattern 
routine was rewritten to accommodate the Yagi-Uda and to 
save storage, and the graphic text segment was changed to 
satisfy the need for additional parameter inputs and outputs. 
The program remains functional for all antennas listed below. 
Figure 21 shows how the parameters and the output appear to 
the operator at the CRT graphics display console. 

To operate the program, first ready the Graphics device 
to be used by following simple lab instructions for setting 
up " Gated . " The program may be quickly loaded into the 
computer by mounting the program " Dump " tape on tape unit #1 
and loading a few cards — BOOT,AJOB, AAGT , AASSIGN X1=MT1A, 
and ARERUN. Next, press the usual operating buttons on the 
9300 console — IDLE, RESET, CLEAR FLAGS, RUN, and CARDS. The 
cards and tape are then read by the computer and within five 
seconds the computer console will print a request for the 
graphics device number that the user desires to interface 
with the computer. The user then responds by typing at the 
IDEV=1* and a carriage return. If a mistake is made in typing 
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start over again. If the message was successfully typed but 
for some reason the text did not appear on the graphics device, 
the quick load procedure may be used by pressing the following 
computer buttons: IDLE, RESET, STEP, RUN; then retype the 

message IDEV=1*. The graphics device should now contain the 
program text, providing all equipment is functional. 

Each non-zero parameter entry must consist of a field- 

length of fewer than four numbers, plus a decimal, e.g. 

parameters should not exceed four characters of information. 

Upon completing a parameter entry, press the graphics TTY 

carriage-return and the blinking-light pointer will sequence 

to the location of the next parameter input. If a parameter 

is zero it need not be typed; press the carriage-return and 

25 

go to the next entry location. Correction to a typing 
mistake may be made only on the four-character field string 
of the parameter where the light-pointer is located. If the 
error is identified after sequencing to the next location the 
program may either be executed as is if no serious program 
error will result, of if the error can cause abort within the 
SDS-9300, either use the fast-restart^ procedure or reload 



2^In addition to the need for a decimal to be entered 
with all data, a blanked-out space (one that has been skipped) 
is the same as a zero. 

^This procedure is similar to the initial start-up: 
type the "Reset Gated" messages at the graphics teletype, and 
at the SDS-9300 press IDLE, RESET, STEP RUN and in response 
to the console message listed by the SDS-9300 type IDEV=1* 
(device #1) . 
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27 

the core "dump" tape again. Either procedure takes less 
than 30 seconds whereas problem solution- time can take longer. 




Figure 21. Photograph of graphic 
display showing parameter input/output. 



27 

This procedure is identical to the initial set-up, in 
that the program dump tape is read into the computer. Press 
the button number "32" at the computer console, and type 
ARERUN at the console. This causes the dump tape to rewind 
and be read in again. Once again the computer should request 
Idev, and the user should respond by typing IDEV=1*. 
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Observe the parameters shown in the photograph. The 
following is a description of these parameters to assist 
the user in specifying his antenna. 

ANTN : Antenna 

1. Tilted Dipole 

2. Vertical Whip 

3. Vertical Whip with Ground Screen 

4 . Inverted L 

5. Sloping Long-wire 

6. Sloping Vee 

7. Horizontal Rhombic 

8. Vertical Half-rhombic 

9. Yagi-Uda 

ELE LENG ; Element length (meters) 

If a single element antenna (ANTN 1 through ANTN 8) 
has been specified, the first available parameter space 
should be used for the element length and the other four 
spaces to the right left blank. If a Yagi-Uda (ANTN 9) has 
been chosen, the following order must be observed: 

ELE LENG 




Directors 

The user may choose to design the array without a reflector, 
in which case the first location at the left must be left 
blank. Any number of directors up to three may be selected, 
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however any blanked directors must be contiguous from the 
right. That is, skipping a director space then entering 
data into the next space, will result in the problem being 
misinterpreted: the first zero director marks the end of 

the array and determines the number of elements present. 

HGT/WDTH : Height and width (meters) 

A single element antenna requires only that the height 
be specified; the last four spaces may be left blank, since 
they will be ignored. If the Yagi-Uda is chosen, select the 
parameters according to the following format: 

HGT/WDTH 

\ 

Height 

The first spacing is for the distance between the reflector 
and driven elements. Subsequent spacings are for the distance 
between respective directors. 



\t 

Spacings 



PHIP : <+> ' Antenna parameter. Normally 4>'=0 will be used 

except under dynamic ship simulation, where the angle will 
specify the array axis relative to the ship's head. See 
figure 5. 

THEP : 0' Antenna parameter. 6'=0° for a vertical YAGI-UDA 

and 0'=9O for horizontal. 

FREQ: Frequency (MHz) 
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EPSL: Dielectric constant of earth 



SGMA : Conductivity of earth (mho/m) 

Typical values are: 



Sea 


1 

i 


Fresh 


! Wet 


. Dry 


Water 


1 

1 


Water 


J Earth 

| 


1 Earth 

1 


80 


1 

1 


80 


1 

| 5-30 


! 2-5 

1 


3-5 


t 

1 

t 


10~ 2 , 10" 3 


1 

! 10 -1 - 10" 3 


1 

0 
1— 1 

0 

r— 1 



SMUL : Multiplier 

1. Multiplies o by .1 

2. Multiplies o by .001 

For numbers < 3 this parameter may be used as a multiplier 
of the diameter of the array elements. The element diameter 
is used in determining the element spacing for computation 
of self impedance. When SMUL > 3 the, expression used for 
spacing is 

d V 2 ” 1 • . 1 i 

, or i. , since r = 1 

V2 200 200. 

(see the computer program instruction that precedes #1925) 

PHI : $ Observation parameter. 

THET : 6 Observation parameter. 

Generally the observer wishes to position himself at 
the point of maximum gain, and take the azimuth and elevation 
sections there. The first solution of the problem will show 
where this point is located by listing the respective angles 
under PK: "Peak" . The user may then enter these values for 

PHI and THET, and solve the problem again. 
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LOG : Log pattern 

If left blank linear patterns will result. If "1." 
is entered, a previously entered pattern via light pen 
will be erased. If "2." is entered, a log pattern will 
result . 

ISTH : Store the horizontal pattern 

ISTV : Store the vertical pattern 

Leave blank except to store the pattern about to be 
calculated for purposes of future comparison with another 
pattern by means of "double-exposure". Enter "1." to 
store. 

IRCL : Recall a stored pattern 

Left blank until patterns have been stored by a 
previous problem solution using ISTH/ISTV. When comparison 
is desired with the solution of the problem being entered, 
enter "1". Ensure that if recall is desired ISTH and ISTV 
are left blank, otherwise comparison will be made of the 
solution with itself. 

ALPH: Array-axis angle 

The tilt angle (or the angle of elevation with the 
ground) of the Yagi-Uda array axis. This is the slope 
angle with the ground for the two sloping antennas (ANTN 5 
and 6) as well as for the vertical half -rhombic (ANTN 8) , 
and the minor angle between elements of the horizontal 
rhombic (ANTN 7) . 
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ISEA: Sea state 



For the tilted dipole, vertical whip, sloping long- 
wire, and Yagi-Uda. Left blank unless iterative solutions 
under dynamic sea conditions are desired. The program is 
written for a sinusoidally varying sea with solutions 
produced at thirty degree intervals. The peak amplitude 
of roll is 8 degrees per sea state, and the peak pitch is 
3 degrees per sea state. 

ICRS : Course of the seas 

Enter the direction of the seas relative to the ship's 
heading . 

GAIN: Gain at the observation points specified by <j> , 0. 

The peak gain, identified under "PK" is also given. 

The <j> and 0 which locate the peak gain are listed adjacent 
to the values specified for the problem (discussed above 
under PHI and THET) . 

FTBR : Front-to-back ratio 

The ratio of the peak forward-azimuth to peak reverse- 
azimuth gains. (This is not a comparison between the 90 
and 270 degree gain values). See footnote 3. 
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REQUIRED FfiR VERT M0N3P6LE WITH GND SCREEN 
FUNCTION A'<EX ( X ) 

COMPLEX akex 
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--ANTN 3 
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ANTN 3 
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• • • fsj y 3 

RQRD FOR V MONO/ V M0N0 WITH SCN> INVERTED L» SL9PING LNG WIRE 
FIJNCTI9N SINC(X) 
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ZZ=(XX««2+K**2*L**2)*#C.5 
ARGlsCMPLX ( 0.0 / -7-Z > 

AR32*C“PLX(0.0/-XX) 
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S4 = SIN(SIGHV-2*K*H*SI\!DL) 
C 0 SUl=Ces<Ul ) 
CCGU2=C8S(U2) 
CeSU3»C0S(U3) 
cesu^t = cos(U4) 
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19A3 C0NJTIMIJE 

V0LDRI =V3LTS( 2) 

IF ( LH ( 1 ) • EQ • 0 « 0 ) V8LTSU ) * V0LDR I ; V8LTS ( 2 ) «0 
CALL SIMO (Z2PA'</V0LTS/2*NE>KS) 

IF (KS.EQ.l) WRITE (6/1950) 
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912 G=(ETHET**2+EPHI**2)/(30*CABS(CURDRI )**2> 

I F ( N • EG • 1 ) GV(I)»G 
IF(N.EQ.P) GH(J)=G 
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REQUIRED FOR YAGI 
SUBROUTINE SIMQ( A, 3 ,N,KS) 
D I YENS I ON A( 1 ),B( 1 ) 
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B( J ) =S A VE/3 IGA 

ELIMINATE NEXT VARIABLE 
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